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£> ' Abstract 

m : 

7— 1 ■ We give explicit expressions for the amplitudes associated with the supersymmetric 

(SUSY) contributions to the process: b — > S7 in the context of SUSY extensions of Standard 
I Model (SM) with non-universal soft SUSY breaking terms. From experimental data, we de- 

O ' duce limits on the squark mass insertions obtained from different contributions (gluinos, 

neutralinos and charginos). 

a: 

Ph! 1 Introduction 

<D . 
43 ■ 

The rare B decays represent a good test for new physics beyond the SM since they are not affected 
appreciably by uncertainties due to long distance effects. Here, in the context of spontaneously 
broken minimal N = 1 supergravity QTJ| , we study penguin diagrams with gluinos, neutralinos 
and charginos, which are responsible AS = 1 radiative mesonic decays. In particular, we study 
the b — > 57 decay [0] that gives the most stringent lower bounds on the average squark mass. 
We know that in generic msugra models [§], the soft universal breaking terms lead to a high 
degeneracy in the sfermionic sector. Flavor changing neutral current (FCNC) tests play an 
important role to constrain the SUSY mass spectrum. 

We thus consider the SUSY extensions of the SM with non-universal soft breaking terms 
We shall use the mass insertion method by which it is possible to obtain a set of upper bounds 
on the off-diagonal terms (A) in the sfermion mass matrices (the mass terms relating sfermion 
of the same electric charge but different flavor) . Obviously the mass insertion method offers the 
major advantage that one does not need the full diagonalization of the sfermion mass matrices. 
Then only a small number of effective parameters (S) summarize the effects. We have applied 
this method to the gluinos, neutralinos and charginos contributions to the decay b —> sj, the 
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charginos contribution being original. From experimental limits, we have then derived upper 
bounds on the off-diagonal terms in the sfermion mass matrices (for squarks down and up). 

This paper is organized as follows. In section 2 we present generalities on the minimal SUSY 
models and the mass insertion method. In section 3, we give the explicit expressions of the 
amplitudes associated with the gluinos, neutralinos and charginos contributions to the decay: 
b — > sj. Finally in section 4, we find explicit expressions for the branching ratio BR{b — > sj) 
and the upper bounds obtained for the off-diagonal terms (A) in the squark mass matrices. In the 
annex, we recall the analytic expressions for the Feynman integrals which arise in the evaluation 
of these amplitudes. 

2 Generalities on the minimal SUSY models and the mass 
insertion method. 

The minimal supersymmetric standard model (MSSM), obtained by supersymmetrizing the SM 
field contents and allowing for all possible soft susy-breaking terms, contains a huge number of 
free parameters. In this note, we concentrate on a specific set of models where these soft-breaking 
terms are close to the msugra universality. 

By minimal supergravity models |J (msugra) we below mean, the low energy limit of spon- 
taneously broken N — 1 supergravity theories which supersymmetrize the SM and present the 
following two features : 

• i?-parity is implemented so that no baryon and/or lepton number violating terms appear 
in the superpotential, 

• the Kalher metric is flat, i.e. all the scalar kinetic terms are canonical. 

These features bring about new sources of FCNC (flavor changing neutral current) effects. The 
experimental limits on B meson physics and in particular b — > 57, constitute interesting FCNC 
tests for these MSS Models, typically requiring squark masses of the same electric charge to be 
relatively degenerate, i.e. their mass difference must be smaller than their average value. Briefly, 
we review the major ingredients which give rise to this new source of FCNC. 
The low energy Lagrangian consists of : 

• The superpotential of the N — 1 globally supersymmetric SM: 



W = h u QH lU c + h d QH 2 d c + h L LH 2 e c + iiE 2 H x (1) 
The scalar part of SUSY soft breaking terms for the minimal N = 1 supergravity theories 



KTft = m 2 x l^l 2 + Am (h u QH lU c + h u QH 2 d c + h L LH 2 e c ^j + Bm[iE 2 E x + h.c 



i=scalar 



(2) 



where A and B are two dimensionless free parameters of the trilinear and bilinear scalar 
contributions; m denote the scale of the low energy SUSY breaking. 
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From equations ([[]) and (|2|), we obtain the squark down 6x6 mass matrix (Q = 



M 2 — 
dd* 



m~ — m~ — 

d L d* L d L d c L 

1 1 

m~ — m~ — 

d* L df did-* 



(3) 



Where 



m i,d* = m kd? = mdTn+d + m2 x 1 ( 4 ) 



and 



m~~ = Amm d + fim d (Hi)/ (H 2 ) (5) 

with m d = 3x3 quarks D mass matrix and Cd = — | is the electric charge. 

At this stage, it is clear that the di — d\ — g coupling cannot lead to flavor changes (FC). 
Indeed, if we diagonalize m d m\ , we diagonalize at the same time m d ~ d n- However, this is no 

longer true if we renormalize m ^ d r'- its value stays ([|) at the high scale, but its evolution down 

to the Mw scale cannot be diagonal due to the h u QH\u c term in the superpotential ([[]). Hence 
the 3x3 mass matrix of rn"-^ renormalized to the M w scale, reads: 

m—^j (g 2 = M 2 ) = m d m d + + m 2 x 1 + cm u m\ (6) 

Where the coefficient c can be computed by solving the set of renormalization group equations for 
the evolution of the SUSY quantities. From equ. (||) we see that the simultaneous diagonalization 
of rn~~ and m d m d + is no longer possible due to the presence of the cm M m+ term. The flavor 

change is proportional to c and to the usual (CKM) angles. In a basis where d L — d\ — lj couplings 
are flavor diagonal, the flavor mixing occurs in the squark propagators. The above remark can 
be summarized in the following schematic way: 



c(V.[mi ia9 } 2 



(7) 



For the L — > R transitions, we have: 



- AL = A^ (8) 



The quantities are mass insertions connecting flavors i and j along a sfermion propagator and 
the indices L, R refer to the helicity of the fermion partners. There are three types of sfermions 



3 



mixing: All, Arr and Alr. In the MSSM case with universal soft SUSY breaking (msugra), 
there exists a kind of hierarchy among mass insertions that is (All)^ 3> (Alr)^ ^> (Arr)^. 
This is no longer true if flavor changing is produced by another kind of "initial" conditions. Then, 
generally, nothing can be said about the hierarchy of these 3 contributions. In that case, one 
needs a model-independent parameterization of the flavor changing (FC) and CP quantities in 
SUSY to test variants of the universal MSSM. The chosen parameterization is the mass insertion 
approximation ||-||. It concerns the most peculiar source of FCNC SUSY contributions that 
do not arise from the mere supersymmetrization of the FCNC in the SM. They originate from 
the FC coupling of gluinos, neutralinos and charginos to fermions and sfermions. One chooses 
a basis for fermions and sfermions states where all couplings of these particles to gauginos are 
flavor diagonal, while the FC originates from non-diagonal sfermion mass terms in propagators. 
Denoting by A the off-diagonal terms in the sfermions mass matrix (i.e. the mass terms relating 
sfermions of the same electric charge, but different flavor), the sfermion propagators can be 
expanded as a series of 5 = A/rh 2 where m is an average sfermions mass and a typical scale of 
the SUSY breaking. As long as the ratio of non-diagonal entries (A) to the average squark mass 
is a small parameter fl5j, the first term in the expansion, obtained from the non-diagonal insertion 
of mass between two diagonal squark propagators, represents an reasonable approximation. This 
method has the advantage that one does not need the full diagonalization of the sfermion mass 
matrices. So, from the FCNC experimental data we may derive upper bounds on the different 
5's. 

In the following section, we give the explicit expression of the amplitudes associated to the 
gluinos, neutralinos and charginos contributions to the decay: b — > 37 pfl , [p70|] , |TT|| . 

3 Amplitudes contributing to the decay b —> 

The MSSM Feynman rules used for the calculation of the amplitudes can be found in the reference 
U . The calculation of these amplitudes is done with these Feynman rules and the mass insertion 
approximation. Supersymmetric penguin diagrams contributing to the decay b — > 57 are: 

• gluinos (pengluinos): figures |l|(a) and |l](b) 

• neutralinos (penneutralinos): figures ||(a) and |2|(b) 

• charginos (pencharginos): figures ||(a) and ||(b); figures £|(a) and |](b); figures |5|(a) and 
|(b). 

These diagrams induce the effective operator Olr = m^e M (g) s(p — q)a^ lJ q u P R b(p), q is the out- 
going momentum of the photon. 

Now we can give the explicit expressions of the amplitudes associated with the supersymmetric 
penguin diagrams. 

3.1 The pengluinos 

From the diagram illustrated in figure |](a), we obtain 
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Figure 1: Gluino contribution (pengluinos). 



T~ = e D C 2 (R)a s ^=V^^e,(q)s(p - q)a^q v P R \pH(X d ) + qH(X~ g )}b(p) (9) 
where Di with i = 1 6 are squark down mass eigenstates, eo is the electric charge of the squark 

M? 

D; Xg = -^r and the function H(Xg) is given in annex; 5ll is the mass insertion connecting 
flavors b and s with the helicity L: 

A (Mh - MpZfrZg All .... 

fal = S K = M (10) 

in which Zd is a mixing matrix defined by: 

*<^> = *S ( Ml R ut ) Z ° ' 

is the average squark down mass and C 2 (R) = | (fundamental representation ).So the T~ 
expression becomes: 
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T~ = e D C 2 (R)a sJ ^=^e,(q)s(p - q)a^ q v P R [pH(X g ) + qH(X~ g )]b(p) 
From the diagram drawn in figure 1(b) we have 

T^ R = e D C 2 (R)a s ^^V^M 1 (X^(q)7(^qJa^q v P R b(j ) ) 

with 

° (Mf, - MpZ^Z^ _ A LR 

LR h M D M D 

Thanks to the experimental limits, it will be possible to put upper bounds on the different 5's, 
that is on the non-diagonal terms in the sfermion mass matrix. 
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HL UiL Q dlt dn 

|(a) 1(b) 
Figure 2: Neutralino contributions (penneutralinos). 





b R U iL U iL ^7 q s L b R 

1(a) |(b) 
Figure 3: Chargino contributions with photon coupling to up squark and mass insertion LL. 



3.2 The penneutralinos 

The penneutralinos are illustrated in figures ||(a) and ||(b). 
• for the diagram (a) 

5ll 



T x° LLj 6Daw 2cos 2 (6 w )M 2 D ^ 



sfa(z Lx °)s{P - q)e^q)a^q u P R [pH(X 0j ) - qH(X 0j )/3]b(p) 



M 2 



j = 1...4 the four neutralinos indices, X 0j = -j-pf with M v o the neutralino mass, and 



^Z 1 ^ sin 9 W - Z% cos 9 V 



Clearly, z Ly a is less or equal to 1. In msugra, for example, we will have: z Lx o 
because Zj^ « 1 for the lightest neutralino ( bino-like), and z Lx o = 0.8 for Z]f « 1 

for the diagram (b) 

sm(9 w )S LR M x o 



sin 2 (9 w ) 



with 



T x° RLj eDaw 3cos 2 (e w )M 2 D ^ 



Z R X ° = G Z N* sin 0™ ~ Z N* cos ^w)Z 



\Za(z Rx o)M 1 (X oj )s(p - q)e^q)a^q v P R b(p) 



2j* 



rlj* 

N 



where, in msugra, we can have z Rx o « 1. 
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1(a) |(b) 
Figure 4: Chargino contributions with photon coupling to the chargino and mass insertion LL. 



3.3 The pencharginos 

Due to the photon-chargino coupling, there are 6 diagrams. For the mass insertion LL, 4 diagrams 
illustrated in figures ||(a) and||(b) and figures |](a-b) contribute. When the helicity flip is realized 
in the quark b external line, only the wino component of chargino is concerned in the calculation of 
the amplitude (diagrams b).But in the case where the helicity flip is realized on the chargino line 
the higgsino components are taken into account (diagrams a). From the penchargino diagrams: 
figures ||(a) and |3](b), where the photon is coupled to the squarks, we obtain: 



for the diagram (a) 



m b 6 LLx .M x7 Z+*Z 2j 



T- = e u a w — *' ' ^M 1 (X j )s(p - q)e fl (q)a^q 1/ P R b{p) 

KLL i M w cos(p)M^v27r 



M 2 



where j = 1,2 are the two chargino states, X,- = and M x - is thej chargino mass, 
e u = 2/3 is the squarks up electric charge, Mjj the squark up average mass. And 

<w, = E E Ea - hK) {M \~ Mh) z$% K zff*v b *j (ii) 



«=1J=1*=1 M u 



in which J and K run over the 3 generations of squarks and Ui with i = 1,...6 are 
up squarks mass eigenstates. As in 5ll, ^LL\j contains squark mixing factors Z , but in 
addition, there are some known Cabbibo quarks mixing factors (e.g. V bc ). 



for the diagram (b) 



T = e u a w * J 2 yfas(p - q)e fi {q)a^q u P R \pH(X j ) + qH(X 3 ) /3]6(p) 

In msugra for the lightest chargino, we have Z±* w 1 when Z^* ~ 0; in such case, we remark 
that only the diagram (b) ( wino component of chargino) contribute to the amplitude. 

The pencharginos illustrated in figures |(a) and [|(b), where the photon is coupled to the chargino, 
give: 



for diagram (a): 



m b 5 LLxj M-Z+*Z;* 



T x - = -Ofc— ; j ' ' ^F(X 3 )s(P - <lK(qWV R b(p) 

Kll j cos(p)M^v27r 
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Figure 5: Chargino contributions with mass insertion LR. 



for the diagram (b) 



$LLxj Z \j Z Xj 



T = a w A ; 2 1J 1J Vas(p - q)s,(q)a^q u P R [pG(X J ) + qG(X,)/2]b(p) 



We define the fraction of gaugino in the chargino j by: z x ij = Z^*Z^. the F(Xj) and 



7 + * , 

G(Xj) are given in annex [| As above, in msugra, only the diagram (b) will contribute to 
the amplitude. 



The LR mass insertion for the chargino contribution are illustrated in figures |](a) and |5|(b). 
Only the higgsino components of the chargino contributes to the amplitude. Therefore, due to 
the sl — xj ~ UiR coupling giving a factor U j Z^ +3 ^Z2j*P r V s j where U J is a Yukawa coupling 
proportional to the associated quark mass, the top quark contribution overwhelms the up and 
charm ones. 

From diagram (a) in figure 5 we thus obtain: 

m t m b S LRXj z Rj M x 
' xZjy " ~ u ~ w M^cos((3) sin(/?)M2 0F 



T v7„. = ~ e uOt w _„,a\^/a\*A r^ V^ M i( X j)s(p - q)s ll (qa ll ' / q I/ P R b(p) 



whereQ 



m 

1=1 u 



and for diagram (b) : 



m t m b S LRx z Rj M 



with: z R j = Z^Z^j*, the fraction of higgsino in the chargino j, its greatest value is 1 and 
the minimum value \ for one of the 2 charginos. In the following section we give the explicit 
expression of the branching ratio for the decay b — > 57 and the upper bounds on the mass 
insertions. 



J It would be desirable to pull out CKM factors from S's, so that they only reflect squark properties. While 
this is arguably possible for Slr x , it requires non-trivial assumptions for Sll x , which is why we kept CKM factors 
in the definition of both. 
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4 BR(b — > 57) Expression and the upper bounds on the mass 
insertions 5. 



The decay b — > 37 is very interested because the rare -B decay represent a good test for new 
physics beyond SM since they are not affected appreciably by uncertainties due to long distance 
effects. 

is 



The branching ratio 1 10 



BR(b -»■ 37) 



BR{B -> X s7 ) r(6 -> s 7 ) 



BR(B -> X c ez7^) r(6 



where 6 — ► cez/ e is dominant then BR(b 57) = r(6 — ► s7)r B .The explicit expression of 
BR(b — > 57) obtained from the calculation exposed in the section || is: 



BR{b -> 57) 



16 7T 2 



m b a s e D C 2 {R) 



Ml 



5llH{X~) 



(13) 



er,a w M Y Q sin z 



9M 2 



0.7 / 



+ 



eY)OL w m h sin 



18 M 2 cos 2 (#J 
a w m b m t M - 
Ml cos(/?) sin^Mt/' 



-5ll#(X, 



0.7 J 



+ 



,F(X,) 



-euM^X^S 



By imposing that each individual term, in the above equation, does not exceed in absolute value 
the experimental data of BR(b — > 57) that is (1 — 4)10~ 4 (which includes the QCD uncertainties 
following flU), we give upper bounds on the different 5. 

We have chosen the values for the supersymmetric particles from the experimental data given 

in 



12[| . Moreover, we have imposed the following conditions 
the average squark mass: M\j 



M D = Mg. 



For the neutralino masses: we choose the LSP mass, with M v o ps 
The lightest stop mass: M~ t > M x o + 30Gev 
The chargino mass: M x - « 2M X? (GUT) 



f- (GUT relation n) 



Otherwise, we lay down Xg 
are defined in the following table : 



M 2 
9 

M 2 



, x 



x" y 



M 2 



Mf 



Y- The others experimental data chosen 



a s 


a 




sm 2 (8 w ) 


BR(b -> 37) 


M w 




m t 




0.12 


1/127.9 


a/4 


0.2315 


1 -> 410" 4 


80.41 


4.5 


170 


1.49 10~ 12 s 



The results are given, in tables |] and [|, for different tan(/3) values ( i.e. tan(/3) =2, 5, 10, 
20 et 40) and for 2 values of M q : M q = 300GeV et Mg = 500GeV (the corresponding results 
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x~ 9 


5ll 


5lr 




X 


5ll z Lx o 




300 


1 


2.96 


io- 2 




50 


3 x 10~ 2 


7.1 


0.34 




(0.36) 


(8.2) 


(2.7 x 10" 






(10~ 2 ) 


(19.7) 


(0.94) 


600 


4 


9.5 


1.8 x 1(T 


-2 


100 


io- 2 


8.4 


0.25 




(1.44) 


(26.4) 


(4.9 x 10" 


- 2 ) 




(4 x 10~ 2 ) 


(23.25) 


(0.7) 


800 


7 


17.6 


2.6 x 10- 


-2 


130 


0.19 


9.8 


0.23 




(2.56) 


(48.8) 


(7.2 x 10- 


- 2 ) 




(7 x 10~ 2 ) 


(27.3) 


(0.64) 



Table 1: Limits on the off-diagonal terms 5ll and Slr for down squarks with Mg = 300GeV (or 
Mq = 500GeV), coming from gluino and neutralino contributions. 



M x - 


X 


&LL X Z X L1 


$lr x z ri 


&LR X Z R1 


$LR x z Rl 


&LR X Z R1 






tan(/5) = 2 


M u=Mq 




tan(/3) = 5 


tan(/3) = 10 


tan(/3) = 20 


tan(/3) = 40 


100 


0.1 


0.57 


0.08 


0.04 


0.02 


0.01 


0.0051 




(4 x 10~ 2 ) 


(1.58) 


(0.09) 


(0.045) 


(0.023) 


(0.011) 


(0.0058) 


200 


0.44 


1.23 


1.6 


0.76 


0.39 


0.2 


0.099 




(0.16) 


(3.41) 


(0.2) 


(0.095) 


(0.049) 


(0.025) 


(0.012) 


300 


1 


2.34 


0.3 


0.14 


0.075 


0.03 


0.019 




(0.36) 


(6.5) 


(0.83) 


(0.4) 


(0.2) 


(0.1) 


(0.052) 



Table 2: Limits on the off-diagonal terms 5ll and 6lr for up squarks with Mg = 300GeV (or 
Mg = 500GeV), coming from charginos contributions. 



are noted between parenthesis in tables) . We have obtained the LL insertion limits by imposing 
BR(b — > sj) = (4) x 10~ 4 and for LR insertion BR(b — > 57) = (2) x 10" 4 to be consistent with 
reference [§] in the gluino case. For chargino contribution we take Mu = Mq. 
From the results obtained we remark that: 

• because of our use of the mass insertions, we are limited to 6 < 1 for the branching ratio 
expression (|13|) to make sense. When larger than one, experimental bounds on 5 like 
5ll < 2.96 thus really mean that the maximal effect of such terms (for 5ll ~ 1) only 
contributes a fraction of about l/9th of the experimental bound; 

• in the gluino case, 5ll is more sensitive to the gluino mass than 5lr, because this last 
contribution has an amplitude enhancement factor of Mg. Then the dependence on M g of 
the H(X g ) function is partially compensated (see the definitions in Annex |^ and the plot 

• for the neutralino, the limit on 5ll is less sensitive to the neutralino mass, thanks to the 
small values of X contained in the same H(X ) function. However, even for the upper 
value for z Lx o = 1 (e.g. for j = 2), the limit cannot be more competitive than the gluino 

limit, except if is smaller than in the msugra models. The 5lr limits decrease weakly 
with the enhancement of M x o, due to the additional power of M x o in the amplitude and 
the fact that the function Mi(Jf ) is fairly constant for small X Q (see Annex || and plot |6|); 
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the chargino contribution is the only one constraining the differences between the up squark 
masses. However the expressions for 5ll x an d 8lr x contain Cabbibo mixing factors. For 
3lr X i the dominating top contribution allows to extract a simple factor of V ts V t b ~ 
The first remark above then applies for interpreting our mass insertions results, once 
> 0.03. Nevertheless, we obtain a limit that quickly becomes more constraining 
with increasing tan(/3). For 5ll X i the factorization of CKM elements cannot be so easily 
justified: if for instance there is large mixing between the squarks 2 and 3, the leading 
CKM factors are of order one, and the expression (|TTD for Sll x with up squarks becomes 



the same as 5ll for down squarks in (|10|) . The sensitivity on Mg of the bounds is about the 
same factor of 3 for the chargino 8ll x as fo r the neutralino 8ll and 5lr while the chargino 
Slr x is less sensitive for small chargino masses; 

the limits in tables (1) and (2) can easily be generalized for different values of squark and 
gaugino masses: they are inversely proportional to the functions given in Annex |6] and 
plotted in figure 

finally, if following |p~3| , we use the latest CLEO data |14]] and substract the SM contribu- 
tion [15|] , we obtain a tighter bound, in the absence of cancellations between the various 
contributions that was assumed throughout this paper. The various results in the tables 
are then reduced by a factor of about 20. 



5 Conclusion 

We have given explicit expressions for the amplitudes associated to the supersymmetric contribu- 
tions to the decay b — > in the context of supersymmetric extensions of SM with non-universal 
soft SUSY breaking terms. The model independent parameterization which we have chosen is 
the mass insertion approximation. From the FCNC experimental data, we have derived upper 
bounds on the different £'s. The contribution from the chargino and neutralino exchanges are 
less sensitive to the gaugino mass than the gluino contribution. 



6 Annex 

-1 + 9X + 9X 2 - 17X 3 + 6X 2 (XlnX + 31nX) 



H{X) 



F(X) 
G(X) 



n(x-r 



5 



Ml(x) = l±*x^*0M±*W = HX)/2 



5 - AX - X 2 + 2 ln(X) + AX ln(X) 
2(X- l) 4 

1 + 9X - 9X 2 - X 3 + 6X(1 + X) ln(X) 
3(X- l) 5 
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0.1 



1 

X 
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Figure 6: The different functions encountered in the evaluation of the amplitudes, for the available 
values of X = M| augino /M s 2 rks on a logarithmic scale. 

Note added 

This work was supported in part by the GdR Supersymetrie of the French CNRS. 

While this work was being refereed, a general study of B — > X s ^ appeared|ll| with interesting 
results on the interferences between various contributions, including the ones presented here. For 
the parameters studied there (/i = 300GeV, M- q = 500GeV, M 2 = 100 -> 230GeV), we agree that 
there is no constrain on up squarks from table [2|: the light chargino is a gaugino with zr\ <C 1 
and the heavy chargino gives the last line in parenthesis, with zr2 = 1- Taking the Cabbibo 
factors into account, the strongest bound on the off-diagonal element S u< lr,23 is just about one 
for tan/? = 50, which means no contraint in the mass insertion logics. 
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